On the Critical Delays of Mobile Networks under 
Levy Walks and Levy Flights 

Kyunghan Lee'^, Member, IEEE, Yoora Kim-1-, Member, IEEE, Song Chong", Member, IEEE, 
Injong Rhee^", Member, IEEE, Yung Yi^, Member, IEEE, and Ness. B. Shroff-f, Fellow, IEEE 
^^{kleeS, rhee}@ncsu.edu, ■'■{kimy, shroff} @ ece.osu.edu, "{songchong, yiyung}@kaist.edu 



O 

CD 



C/2 



X: 



Abstract — Delay-capacity tradeoffs for mobile networks have 
been analyzed through a number of research work. However, 
Levy mobility known to closely capture human movement pat- 
terns has not been adopted in such work. Understanding the 
delay-capacity tradeoff for a network with Levy mobility can 
provide important insights into understanding the performance 
of real mobile networks governed by human mobility. This paper 
analytically derives an important point in the delay-capacity 
tradeoff for Levy mobility, known as the critical delay. The 
critical delay is the minimum delay required to achieve greater 
throughput than what conventional static networks can possibly 
achieve (i.e., 0{l/y/n) per node in a network with n nodes). The 
Levy mobility includes Levy flight and Levy walk whose step 
size distributions parametrized by a G (0, 2] are both heavy- 
tailed while their times taken for the same step size are different. 
Our proposed technique involves (i) analyzing the joint spatio- 
temporal probability density function of a time-varying location 
of a node for Levy flight and (ii) characterizing an embedded 
Markov process in Levy walk which is a semi-Markov process. 
The results indicate that in Levy walk, there is a phase transition 
such that for a £ (0, 1]> the critical delay is always 0(n2 ) and 
for a G [1, 2] it is 0(n"2 ). In contrast. Levy flight has the critical 
delay e(nt) for a £ (0,2]. 



I. Introduction 

Since the seminal work by Gupta and Kumar [IJ on the 
capacity of wireless networks, delay and throughput tradeoffs 
for wireless networks have been extensively studied for vari- 
ous mathematical techniques, scheduling algorithms, channel 
models, mobility models and physical layer techniques. The 
work by Grossglauser and Tse f2| showed that the per-node 
throughput remains constant (0(1)) when node mobility is 
used for communication. This result is surprising because 
Gupta and Kumar IT] had previously shown that the per- 
node throughput {0{1/ y/n)) in wireless networks with no 
mobility diminishes as the number of nodes n increases. This 
throughput gain is achieved at the cost of larger delays. 

The amount of delay that a network needs to sacrifice to 
guarantee a given throughput has been studied under various 
mobility models ll3l-||5]. In particular, Sharma et at. |j6] 
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(a) Brownian motion (b) Levy mobility (c) Random waypoint 

Fig. 1. Sample trajectories of (a) Brownian motion, (b) Levy mobility and 
(c) random waypoint. 



Studied the minimum delays required to achieve more per- 
node throughput than 0(l/A/nlll under various mobility mod- 
els including i.i.d., random waypoint, random direction and 
Brownian motion. This minimum delay is called critical delay. 
However, although the work provides a nice framework for 
studying delay-capacity scaling for wireless networks under 
a family of random walk models, the practical values of 
these mobility models are limited. While these models are 
simple enough for mathematical tractability, they do not reflect 
realistic mobility patterns commonly exhibited in real mobile 
networks. 

Humans are a major factor in mobile networks as most 
mobile nodes or devices (smartphones and cars) are carried or 
driven by humans. Recent studies Q-lHl on human mobility 
show that step size distribution^ are heavy-tailed where a step 
is defined to be the straight line trip of a moving object (e.g., 
particles or humans) from one location to another without 
a directional change or pause. These mobility patterns are 
accurately modeled by Levy process fTU\. 

Levy mobility is a random walk mobility whose step size 
distribution is parametrized by a G (0, 2] and is heavy-tailed 
except in the extreme case of a = 2. For a G (0, 2), the 
distribution is well approximated by a power-law distribution 
l/z^"*"" where z is a step size. For a = 2, the step size 
conforms to Gaussian distribution!! Intuitively, such a random 
walk contains many short steps and a small yet significant 
number of exceptionally long steps. With different values of a, 
the movement patterns of Levy mobility models are widely 
different. Smaller a induces a larger number of long steps. 

'As Q] showed, @{l/y^) is the maximum throughput that wireless 
networks relying on naive multi-hop transmissions can achieve without the 
help of node mobility. 

-Step size is often referred to as flight length in some literatures. 

^Levy mobility becomes Brownian motion in the extreme case of « = 2. 
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This type of mobility patterns is significantly different from 
Brownian motion and random waypoint as illustrated in Fig.[T| 
In the literature, there are two types of Levy mobility models 
for classification: Levy flight and Levy walk. In Levy flight, 
every step takes a constant time irrespective of its step size 
and in Levy walk, it takes a constant velocity. Levy flight and 
Levy walk can show the same pattern of traces but their time 
durations taken to have such traces are essentially different. 
Intuitively, Levy flight can be easily slotted while Levy walk 
is not. 

Unfortunately, understanding tradeoffs between throughput 
and delay under Levy mobility models is technically very 
challenging and underexplored. Unlike the other random walk 
models permitting mathematical tractability, the Levy process 
is not very well understood mathematically despite significant 
studies on Levy process in mathematics and physics. Thus, 
the conventional techniques JSj, 161 used to study delay- 
capacity tradeoffs cannot be applied to Levy mobility mod- 
els, especially to Levy walk which has high spatio-temporal 
correlation. In more specific, since Levy walk is not eligible 
for discretization for Markovian analysis, its mathematical 
characteristics such as joint spatio-temporal probability density 
function (PDF) are hardly known. Due to such a difficulty, 
analyzing Levy walk is generally considered to be very 
challenging. 

Our main contribution is to analytically derive important 
tradeoffs between delay and capacity for both Levy mobility 
models. An important point in this tradeoff is the "critical 
delay" which is the minimum delay for a mobile network 
to obtain a larger throughput than Q{l/y/n). Our technique 
involves (i) analyzing the joint spatio-temporal PDF of a time- 
varying location of a node and the diffusion equation of the 
node for Levy flight and (ii) characterizing an embedded 
Markov process inherent in Levy walk which is a semi-Markov 
process. Since a different value of a induces a different 
mobility pattern, it also induces a different critical delay. 
Below we summarize our main results. 
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Given that many human mobility traces are shown to have 
values of a between 0.53 and 1.81 Q, according to our results, 
mobile networks assisted by human mobility have critical 
delays between 0(n''-^'') and 8(?7.°'^^). Note that our results 
give much more detailed prediction of the critical delay for 
such mobile networks depending on a while Brownian motion 
and random waypoint always show 0(n) and Q{n'~'-^) for their 
critical delays 161. 

The rest of the paper is organized as follows. We first 
overview a list of related work in Section and introduce our 
system model in Section |III] More details of Levy mobility 
model parameterized by a are described in Section HV] and 
the critical delays under Levy flight and Levy walk are 
investigated in Sections |V] and IVT] respectively. Finally, we 
provide a high level interpretation of our main results in 
Section IVIII and concluding remarks in Section IVIIII 



II. Related Work 

Gupta and Kumar H] showed that the per-node throughput 
of random wireless networks with n static nodes scales as 
a function of 0(1 /^/n) and proposed a scheme achieving 
Q{l/y/n \ogn). The result for static wireless networks was 
later enhanced to 9(l/yn) by exercising individual power 
control im . lfT2l . Grossglauser and Tse l2| proved that a 
constant per-node throughput is achievable by using mobility 
when the nodes follow ergodic and stationary mobility mod- 
els. This result disproved the conventional belief that node 
mobility can negatively impact network capacity as it causes 
connectivity breakup and channel quality degradation. 

Many follow-up studies lH, 0, ini-lnl have been de- 
voted to understand, characterize and exploit the tradeoffs 
between throughput and delay. Especially, the delay required 
to obtain the constant throughput 0(1) has been later studied 
under various mobility models IH, lT6l - llT8l . The studies 
provided that the delay to obtain 0(1) of per-node throughput 
becomes Q{n) for most mobiUty models such as i.i.d. mobility, 
random direction, random waypoint and Brownian motion 
models. 

Another interesting question that has attracted researchers is 
what should be the minimum delay to achieve asymptotically 
higher throughput than 0(l/\/n), the per-node throughput of 
static networks. This has been studied under the notion of 
critical delay ||5], for two families of random mobility 
models: hybrid random walk and random direction. The hybrid 
random walk model spUts the network of size 1 with n^^ cells 
and further splits a cell into n^~'^^ subcells for /? G [0, 1/2]. 
Then, a node moves to a random subcell of an adjacent cell in 
every unit time slot. In this model, i.i.d. mobility corresponds 
to = and random walk mobility corresponds to /3 = 1/2. 
For any [5 G [0, 1/2], critical delay is proved to be Q[v?^). The 
random direction model chooses a random direction within 
[0, 2tt\ and moves to the selected direction with a distance of 
'nT"' with a velocity n^^/^ for 7 G [0,1/2]. In this model, 
random waypoint and Brownian motion are represented with 
7 = and 7 = 1/2, respectively. The critical delay is proved 
to be 9(^1/2+7). 

III. Model Description 
A. System Model 

We consider a wireless mobile network indexed by n, 
where in the n-th network, n nodes are distributed uniformly 
on a completely wrapped-around square S{n) whose width 
and height scale as ^/n and the density is fixed to 1 with 
increasing ?70 Without loss of generality, we set the width 
and the height of the square S{n) as ^Jn. We assume that 
all nodes are homogeneous in that each node generates data 
with the same intensity to a per-source destination. The packet 
generation process at each node is assumed to be independent 
of node mobility. 

A source-to-destination packet can be delivered by either 
direct one-hop transmission or over multiple hops, say k hops, 

^This model is often referred to as an extended network model. In another 
model, called a unit network model, the network area is fixed to 1 and the 
density increases as n while the spacing and velocity of nodes scale as 1/ y/n. 
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using relay nodes. We call it fc-hop relay transmission. We 
assume that all nodes can serve as relay nodes for other source 
nodes and the nodes serving as relay nodes can only forward 
packets rather than replicating packets (not to overproduce the 
same packets in the network). 

To model interference in wireless networks, we use the 
protocol model as in |[T|, under which nodes transmit packet 
successfully at a constant rate W bits/sec, if and only if the 
following is met: let Xi{t) (g M^) denote the location of 
node i{i = l,...,n) at time t{> 0). For a transmitter i, 
a receiver j and every other node k 7^ i,j transmitting 
simultaneously, 

d{Xk{t),X,{t)) > {l + A)d{X,{t),X,{t)), 

where d{x,y) denotes the Euclidean distance between loca- 
tions x,y E R^, and A is some positive number 

A packet can be delivered through a scheduling scheme 
which consists of replication or forwarding. We assume that 
only source nodes replicate packets and all other relay nodes 
forward them. As the names imply, replication copies a 
packet and the packet transmitter keeps the packet, whereas in 
forwarding the transmitter does not keep the original packet 
after successful transmission. This selective replication and 
forwarding depending on the node type are often applied to 
suppress the overflow of redundant packets in the network. 
Packets are delivered in two ways: neighbor capture and multi- 
hop capture. In neighbor capture, using mobility, relay or 
source nodes are located within the communication range of 
the destination. In the multi-hop capture, a source establishes a 
multi-hop path to the destination and delivers the packets over 
the path. We assume a fluid packet model ||T9l so that the 
delivery can occur immediately even in the case of multi-hop 
capture because the transmission delay is negligible compared 
to the delay from node mobility. We denote by 11 the class of 
all scheduling schemes conforming the descriptions above. 

B. Performance Metrics 

The primary performance metric in many networking sys- 
tems is per-node throughput measured by the long-term aver- 
age of received packets aggregated over nodes: 

Definition 1 (Per-node throughput): Let A7r(n) denote the 
per-node throughput in the n-th network under a scheduling 
scheme tt e 11. It is then given by 

Att ('T-) — lim inf — ^''^ , 

1=1 

where Xyr-.tit) is the total number of bits received at a desti- 
nation node i up to time t under 7r|f| 

Another important metric is average delay: 

Definition 2 (Average delay): Let Dt^^ti) denote the aver- 
age delay in the n-th network under a scheduling scheme 
TT € n. It is then given by 

^ n ^ k 

i=i j=i 

^For simplicity, we omit the subsciipt vr in A^(n) unless confusion arises. 



where is the individual packet delay that a packet j 

experiences to arrive at a destination node i from its source 
node under tt. 

We give special attention to the notion of critical delay, first 
introduced in ||6l: 

Definition 3 (Critical delay): The critical delay in the n- 
th network, denoted by Cn(n), is the minimum average delay 
that must be tolerated under a given mobility model to achieve 
a per-node throughput of uj{l/y/n), i.e., 

Cn{n) = inf £',r(?^)■ 

{7ren:A,(Ti)=w(l/v^)} 

Per-node throughput 6(l/\/?i) is achievable by a schedul- 
ing scheme in static multi-hop networks U. Since node 
mobility can increase per-node throughput at the cost of larger 
delay, the critical delay quantifies the amount of delay that a 
network should sacrifice to achieve the guaranteed "baseline" 
per-node throughput. It can be used as a simple, yet useful 
metric for a mobility model, representing how sensitive the 
delay is to increase per-node throughput. 

Computing critical delay consists of multiple steps. We 
start by following the initial step in ||5l, ||6l which connects 
critical delay to the first exit time. Let denote a disc 

within the square S{n) whose radius scales as Q{y/n). Critical 
delay can simply be regarded as the maximum time duration 
that a node cannot exit from the disc 'D{n) with probability 
approaching 1 as n goes to 00. In our extended network 
model, the average distance from a source node to a destination 
node is <d{y/n) when they are uniformly distributed on S{n). 
Therefore, if nodes travel up to a distance 0{^/n), for a certain 
time duration, the distance between a source or a relay and 
a destination still remains on average which results 

in 0(1/ y/n) per-node throughput (see Lemma [Til. Thus, it 
is obvious that a network aiming at obtaining uj{l/^/n) per- 
node throughput must allow a delay which is no less than the 
maximum time duration that the first exit of a node from the 
disc does not occur with probability approaching 1. This 
insight can be formally described with the notion of the first 
exit time: 

Definition 4 (First exit time): Let Xi{0) = x. The first exit 
time for a disc of a radius r, denoted by T(r), is defined as 

Tir)^M{t>0:X,{t)(^Bix,r)}, 

where B{x,r) denotes the set of points y in S{n) such that 
d{x,y) < r. 

Without loss of generality, we set the radius of the disc 

as Cd\fn where is a constant in the range (0, 1/2). Then, 

critical delay Cx\(n) can be obtained by 

Cn(n) = sup |i(ri) : lim P{T(cdVn) > ^W} = l| • 

Lemma 1 (^, ^): Suppose that on average each packet 
is relayed over a total distance no less than Q{^/n) in an 
extended network model. Then, the per-node throughput A(n) 
scales as 0(l/\/n). 
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IV. Mobility Models: Levy Flight and Levy Walk we have 



In this section, we formally define Levy mobility model: 
Levy flight and Levy walk. 

Levy flight and Levy walk processes are treated separately 
in the literature Il20l - ll22l . Levy flight takes a constant time for 
any step irrespective of its step size, whereas Levy walk takes 
a constant velocity for every step. Thus, in Levy walk, the 
time taken for each step is proportional to the step size. The 
distinction between Levy flight and Levy walk is often made 
based on the speeds of their actual processes. Levy flight is a 
"fast" mobility model which can reach its next destination in a 
constant time no matter how far it is. In a similar context. Levy 
walk falls under a "slow" mobility model. An experimental 
velocity model suggested as a function of step size in fT] 
verifies that a human mobility lies in between Levy flight and 
Levy walk. For convenience, we use Levy mobility model to 
indicate both of Levy flight and Levy walk, unless explicitly 
stated. 

Let Z be a random variable denoting the step size under 
Levy mobility model. Then, Z is generated from a random 
variable Z having the Levy a-stable distribution 1231 by the 
relation Z = \Z\. The PDF of Z is give by 



1 r 



(1) 



where if^it.) = E[e* ] is the characteristic function of Z 



and is given by ip^i^) 



Here, IcI > is a scale 



factor which is a measure of the width of the distribution, 
and a G (0, 2] is a distribution parameter and specifies the 
shape (i.e., heavytail-ness) of the distribution. The step size Z 
for a G (0, 1) has infinite mean and variance, while Z for 
a G [1,2) has finite mean but infinite variance. For a = 2, the 
Levy a-stable distribution reduces to a Gaussian distribution 
with zero mean and variance cr^ = 2c^, and consequently the 
step size Z has finite mean and variance. 

Due to the complex form of the distribution, the Levy a- 
stable distribution for a e (0, 2) is often given as a power-law 
type of asymptotic form, closely approximating the tail part 
of the distribution 123)1: 



1 



l+a ■ 



(2) 



For mathematical tractability, in our analysis we use the 
asymptotic form ^ instead of the exact form ([T]i for a E 
(0,2) while using the exact form ([T]i for a = 2. The form 
(|2]l is known to closely approximate ([T]i and several papers 
in mathematics and physics, e.g., Il20l . l24l . analyze Levy 
mobility using form (|2]i. For the range of Z, since we use 
the extended network model, the step size Z is assumed to 
have a lower bound at 1 and an upper bound at ^/n, i.e., 
P{1 < Z < = 10 Thus, the complementary cumulative 
distribution function (CCDF) of Z becomes P{Z > z} = 1 
for z < 1 and P{Z > z} = for z > ^/n. For z e [1, ^/n), 

'^The bounds ai'e chosen equivalently to the lower bound at 1/y/n and the 
upper bound at 1 for the step size in the unit network model (6). 



P{Z>z} = 



(v^) 



for a e (0,2), 



(3) 



Here, erf(-) is the error function defined as erf(a;) = 
exp(— t^) dt, and c{n) is defined a^ 



c(n) = 




for a e (0,2), 
for a = 2. 



Note that as n goes to oo, the CCDF P{Z > z} for z > 1 
goes to 1/z" for a e (0,2). 

In our analysis, we use the following assumptions on the 
Levy mobility model: (Al) the time taken for each step in 
the Levy flight is set to 1, and (A2) the velocity taken for 
each step in Levy walk is set to 1. Note that as long as these 
two metrics are constant, the scaling property of critical delay 
remains the same, which justifies our assumptions. 

V. Critical Delay Analysis for Levy Flight 

In this section, we will show that the critical delay C^(?^) 
under Levy flight with a distribution parameter a <S (0,2] 
scales as 8(71^) (Theorem [T]). In Section [V-AI we explain 
technical challenges and our approach for proving Theorem [T] 
In Section lV^ we prove Theorem[T]by showing that the upper 
bound on Cn{n) scales as 0(n"2) (Lemma |3]l, and that the 
lower bound on Cn{n) scales as ^l{n~) (Lemma |4]i. 



A. Technical Approach 

We begin with deriving a relation between the first exit 
time of a 2-dimensional random process and the one for its 
1 -dimensional projected process. We then describe trapping 
phenomenon in a diffusion process that have a direct connec- 
tion to the first exit time of a 1 -dimensional random process. 

It is clear from Definition |4] that the statistical proprieties 
of the first exit time do not depend on the choice of node 
index i. Thus, we omit the node index i in the rest of the paper 
Denote X{t) = {Xx{t), Xy{t)) and consider the projected 
processes {Xx{t)}t>o and {Xy{t)}t>o onto i-axis and y-axis, 
respectively. We define for the projected processes the first exit 
time similarly to Definition |4| 

Tx{r) ^ inf {t > : |X,(i) - X,(0)| > r} , 
Ty{r)^M{t>0:\Xyit)-XyiO)\>r}. 

Since the event {|Xa;(<) — ^^^(O)! > r} implies the event 
{d{X{t),X{0)) > r}, we obtain 



P{r.(r) <t}< P{T{r) < t}. 



(4) 



'To be precise, c(n) is also a function of a, i.e., c(n) = c(n, «). Since 
we focus on scaling properties with respect to n for a fixed a, we omit the 
argument a in c(n, a) for notational simplicity. By the same reason, in the 
rest of the paper, we emphasize only n in all variables that depend on both 
n and a. 
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In addition, it is clear that 

P{T(r) <t}< P{T,{r/V2) < t or Ty{r/V2) < t} 

<2P{T4r/V2) <t}, (5) 



where the second inequality comes from the union bound and 
the symmetry of node motion. Combining ^ and (|5]), we have 
for all t > 0, 



HT.ir) <t}< P{T(r) < t} < 2P{T,{r/V2) < t}. (6) 

Our technical approach is mainly based on (|6]l, and is to 
bound the first exit time distribution for 2-dimensional Levy 
flight by the one for the corresponding 1 -dimensional projected 
process {Xx{t)}t>o. We henceforth study the first exit time 
distribution for the process {Xx{t)}t>o- 

The first exit time analysis for 1 -dimensional random pro- 
cesses has been intensively studied in physics and mathemat- 
ics, e.g., Il25l . Specifically, trapping phenomenon (of a diffus- 
ing particle) in physics and its related theories have a direct 
connection to our first exit time problem as explained in the 
following; consider a particle that diffuses in a finite interval 
[0, 2r] (c M) having trapping boundaries at x = 0, 2r. Let 
L{t) (s M) be a random variable denoting the location of the 
particle at time t. The particle is assumed to be initially located 
at L(0) = r, and eventually it is trapped at either of both 
boundaries with probability I. Upon the particle is trapped, 
it disappears in the interval. We call the state of the particle 
survival state until the particle is trapped and disappears. By 
convention, we let L{t) = if the particle is not in survival 
state at time t. If we assume Xx{0) = L{0) r), then Xx{t) 
and L{t) for t > are related as follows: 



Lit) 



d^lX^it) if i<T,(r), 
|0 ift>T,ir), 



(7) 



where = denotes "equal in distribution" 
from dTji that 



Hence, we have 



P{T,{r) <t} ^P{L{t)^0}. 



(8) 



That is, the survival time of a particle in the trapping model 
has the same distribution as the first exit time Tx{r) of a node 
under Levy flight. 

The technical approach for analyzing the critical delay in the 
literature is as follows. In the case of Brownian motion, there 
are two general techniques in studying the critical delay. One is 
to discretize mobility and then apply a Markovian analysis fjj. 
The other is to use a continuous mobility model and solve a 
diffusion equation to obtain a joint spatio-temporal PDF of 
a time-varying location of a node [Sj. The latter enables one 
to obtain the distribution of L{t) whose spatial derivative is 
often referred to as occupation probability^ The occupation 
probability of Brownian motion can be decomposed to find the 
components constituting it. From this decomposition process, 
we find that there is a dominating term which characterizes 
the limiting behavior of the first exit time distribution. 

*The occupation probability in a trapping mode! corresponds to the joint 
spatio-temporal PDF in a random walk model. The mathematical definition 
and the distinction between the occupation probability and the joint spatio- 
temporal PDF will be given in Section IV-BI 



In the case of Levy flight, the joint spatio-temporal PDF 
has a similar form to that of Brownian motion. In addition, 
the occupation probabilities and the first exit time distributions 
for Brownian motion and Levy flight have similar structures 
in the aspect of the dominating terms. Hence, by identifying 
and characterizing the dominating term for Levy flight, we can 
obtain the critical delay under Levy flight. 

B. Analysis 

In this subsection, we provide the detailed result for the 
critical delay under Levy flight. Our main result is derived by 
following three steps: (i) the occupation probability is obtained 
from the solution of a differential equation that governs the 
movement of a particle, (ii) From the occupation probability, 
we obtain the survival probability (which will be defined 
later), which in turn yields the first exit time distribution, (iii) 
By investigating the limiting behavior of the first exit time 
distribution, we can finally obtain the order of the critical 
delay. 

Step 1: Let P{x,t) = j;P{L{t) < x}. Intuitively, P{x,t) 
represents probability that the particle is located at x at time t. 
We call P{x, t) the occupation probability, and it has the 
following properties: 

. (PI) limf^oo Pix, t) = V.T e M. 

. (P2) J^'' P{x, 0) d.T = P{L{0) = r} = 1. 

. (P3) /o P{x, t) dx < 1 Vt > 0. 

. (P4) P(0, t) = P{2r, = Vt > 0. 

• (P5) Since P{x, 0) is a PDF having a support {r}, we 
have P{x,0) = Sx,r, where Sxi,x2 denotes the Kronecker 
delta which is defined to be 1 if a;i = X2 and otherwise. 

To be precise, P{x,t) for t > could not be a PDF due to 
(P3). However, the function obtained by normalizing P{x,t) 
with the integral J^^ P{x, t) dx, denoted by P{x, t), becomes 
a PDF for a finite time t. We call P{x, t) the joint spatio- 
temporal PDF at location x and time t. 

In the first step, we obtain the occupation probability P{x,t) 
for the process {L{t)}t>o. For this, we need to characterize 
the associated 1-dimensional process {Xx{t)}t>o- We first 
consider the case of a € (0, 2) and summarize the result in 
the following lemma. 

Lemma 2: Suppose that {X{t)}t>Q is 2-dimensional Levy 
flight with a distribution parameter a G (0, 2). Then, as n 
goes to oo, the projected process onto x-axis {Xx{t)}t>o 
approaches to 1-dimensional Levy flight having the same 
distribution parameter a. It holds for the process {Xy{t)}t>o. 

Proof: Let Zi and 6i {i — 1,2, . . .) be random variables denot- 
ing the i-th step size and direction of the process {X{t)}t>o, 
respectively. Then, X{t) for t = 1, 2, . . . can be expressed as 



Xit)^iXx{t),Xyit)) 



x(o) + 



(9) 



We will show that, as n goes to oo, arbitrary step size of 
the projected processes (i.e., Zi\cos6i\ and Z^j sin 6*^1) has a 
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power-law type CCDF with an exponent a, i.e., for z > 1, 

lim P{Z,| cos 6*, I > z} = lim P{Z,\ sine,\ > z} 

71— >-oo n-^oo 

_ c* (10) 



where c* = -I J^^ (cosi?)"d?9. Since the projected processes 
take a constant time for every step irrespective of step size, 
the property in ( fTOl i proves the lemma. 

Now we prove (fTOl i. By conditioning on the values of the 
random variable 6*; ~ Uniform[0, 27r], we can rewrite the 
CCDF of Zi\ cos 0i\ as 

P{Z,| cose,\ >z}^ P{Z,\ cose,\ >z\0,^ ^}dFeM 
Jo 
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1 

2 



P{Zj|cosz?| > z}di9 



P{ZjCosi? > z}di?, 



(11) 



where the last two equalities come from the independence 
of the random variables Zi and Oi, and the symmetry of 
the function |cosi9|, respectively. Using (O, the probability 
V{Zi cos I? > z} in ((TT) can be obtained for d G [0, |] as 

P{ZiCosz? > z} 

^ f c{n) ■ {{^^r ~ {^T) for ^ e [O,cos-i(^)), 
Hence, the CCDF P{Zi| cos 6*^1 > z} is given by 



P{Z,|cos0,| > z\ 



(cost9)"d7? 



2c(ri) 



(12) 



Noting lim„^oo c{n) = 1 and lim„ 
have from ( fT2] i that 



' (^) 



K we 



lim P{Z,| cos6lj| > z} — / ' (cosi9)"di? = — . 



Since | siii(?i| = | cos 6*^1 for ^ Uniform[0, 27r], we have 
P{Z,|sin(?,| > z}-P{Z,|cos^^,| > z}, 

which completes the proof. □ 

Motivated by Lemma |2] and (|7]i, we study the occupation 
probabiUty for 1-dimensional Levy flight with a G (0,2) 
in a finite interval [0, 27-] having trapping boundaries. For 
mathematical tractability, our study in this subsection assumes 
continuous limit where the scale factor |c| in ([TJ approaches to 
zero. Then, the occupation probability P{x,t) for a £ (0, 2) 
is governed by the following fractional Fokker-Planck equa- 
tion Hal Eq. (22)], im Eq. (28)]: 



dP{x,t) d°'P{x,t) 



dt 



9|a;|" 



(13) 



where F {= > 0) is a generalized diffusion coefficient and 
is the Riesz-Feller derivative of fractional order a 1271 . 



scale factor |c| approaches to zero, the 2-dimensional Levy 
flight converges to a Wiener process which mathematically 
models a continuous movement of Brownian motion. Since 
1-dimensional projected process of 2-dimensional Brownian 
motion is also Brownian motion ||5l, the occupation probability 
for a = 2 is governed by the normal diffusion equation where 
the spatial derivative of order a with a G (0, 2) in ( fTsT i is 
replaced by the second order derivative with a = 2 ||25]| . 
Therefore, with continuous limit, the occupation probability 
P{x,t) for a G (0,2] can be described by the differential 
equation in (fT3T l. Through Appendix |A] we show that the order 
of the critical delay under Levy flight does not change with 
continuous limit. 

Applying the standard method of separation of variables 
gives the solution of (fTsl i as follows: 



P{x, t) = Y^ h.ip, (x) cxp {X,Ft) . 



(14) 



1=1 



Here, hi {i ~ 1,2, ...) aie determined from the initial condi- 



tion P{x, 0) 



(as shown in (P5)) and are given by hi 



ipi{r). The functions ipiix) and the constants can be ob- 
tained from the solutions of the problem D[?/'i(.T)] = Xiipi{x) 
for the operator T) = j^j^' ^nd are called eigenfunctions and 
eigenvalues of S, respectively. Without loss of generality, we 
assume that are arranged as |Ai| < IA2I < • • • . 

Step 2: Let S(t) = P{L(t) 7^ 0}. Intuitively, S{t) 
represents probability that the particle has not hit any trapping 
boundary by time t. We call S(t) the survival probability. 
The survival probability can be obtained from the occupation 
probabiHty P{x, t) by S{t) = /J'' P[x, t) Ax. Thus, from (O, 
the survival probability is given by 

S{t)^y2i',{r) tl;,{x)dxexp{KFt). (15) 
i=i Jo 

The first exit time distribution P{Tj.(r) < t} can be obtained 
from the survival probability S{t) through the following 
relation: 



nTx{r) <t} = P{L{t) = 0} = 1 - S{t). 



(16) 



Here, the first equality comes from ([8]) and the second equality 
comes from the definition of S{t). By combining ( fTSl ) and 
( fTSl l, we obtain the first exit time distribution in terms of the 
eigenfunctions ■4ji{x) and the eigenvalues Ai as follows: 

00 „2r 

HT.{r) <i] = l- VV'^(0 / i^^{x)dxe^p {\Ft) . (17) 

^=l J" 

For a = 2, the eigenfunctions and the eigenvalues in ( fTTb 
can be obtained from the boundary conditions P{0,t) ~ 
P{2r,t) = Vt > (as shown in (P4)), and are given by 

- \/^sin(fr) ™d A, = - (f respectively El. 
For a G (0,2), Gitterman |26l provided a solution of ( fT3] ) 
whose eigenfunctions and eigenvalues are given by V'«(a;) = 
^"'n(^) and A^ = ^ respectively. Thus, under 



d\x\ 

We next consider the case of a = 2. In this case, as the Levy flight with a G (0, 2], the first exit time distribution can 
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be expressed as an infinite series of exponential functions as 
follows: 



oo 



1 



(18) 



where rji 



J\ 2{1 — cos(i7r)} 



sin(f) and ^ F(f ) 



As will be shown later in the proof of Lemmas |3] and |4] 
the smallest (i.e., dominant) decay rate in the exponential 
functions in ( fT4l i (i.e., | Ai |) determines the limiting behavior of 
the first exit time distribution. That is, the smallest decay rate 
characterizes the critical delay under Levy flight. The solutions 
in ||251 . Il26l show that the dominant decay rate |Ai| scales as 
e(r-") for a G (0,2]. 

Step 3: We are now ready to derive the main result of 
this subsection. By using the closed-form expression for 
P{Tx{r) < t} in ( fTsl l, we investigate the order of the critical 
delay, stated in Lemmas |3] and |4] 

Lemma 3 (Upper bound for Levy flight): Suppose that un- 
der Levy flight with a distribution parameter a S (0,2], the 
time t = i{n) in ^{T{cd\/n) > t} scales as 6(n7+'^) for an 
arbitrary e > 0. Then, we have 

lim P{T{cdV^) > i{n)} = 0, 

n— >oo 

which shows that the critical delay Cn{n) under Levy flight 
scales as 0(nt). 

Proof: We will prove this lemma by showing that 
liuin^co P{Tx{cd\/n) < i{n)} = 1. Then, by substituting 
r = Cd\/n and t = i{n) into ^ and taking a hmit to n, 
we obtain 

1 = lim P{Tx{cdV^) < i{n)} 

n— ^oo 

< lim P{T{cdVn) < i{n)}. 

n— f C30 

That is, \iinn^ooP{T{cd\/ri) < t{n)} = 1, or equivalently, 
liuin^oo P{T{cd^/n) > t{n)} — 0, which proves the lemma. 

First, consider the case of a = 2. We substitute r — Cd\fn 
and t = t(n) into ( fTSl l. Then, the series on the right-hand 
side of ( fTsT i becomes a function of n, and (for notational 
convenience) we let 



i(n) 



P{r,(cdV^) < On)} = 1 - ^ cxp ( - 

4=1 

= l-5(n). 



{cd)'^n 



We now need to take a limit to S{n). To validate the 
interchange of the order of limit and summation, we will 
show that there exists a constant h E N such that the 
infinite series S{n) converges uniformly onV = [n, c>D)0The 
uniform convergence will be shown by using the well-known 
Weierstrass M test 1281 . 

Since t{n) — 8(n^+'^), there exist constants h G N and 
c > such that 



t(n) > cn 



l+c 



for all n > h. 



(19) 



Let TO = F{T:/2cdfc{ny (> 0). Then, flie i-\h function of flie 
series Sin) is bounded by a constant Mi = ^{exp(— m)}' for 
all n > fi as follows: 



r]t exp 



{cdYn 



tin) 



4 

< — exp 
4 

< — exp 

TT 



icd)^n 



,l + £ 



c(n)< 



Here, the first inequality comes from the bounds |?7i| < 
^ Vi G N and ( fT9] l, and the second inequality comes from 
the bounds P > i yi & N and n"^ > (n)'^ Vn > fi. Note that 
the series J^'iLi ^'^i converges since it is a geometric series 
with a common ratio cxp(— to) G (0,1). Since the target of 
the functions is a complete normed vector space, the infinite 
series S{n) converges uniformly on V. Consequently, we can 
interchange the order of hmit and summation, and we have 

lim P{T,(cd^^) < t(n)} 

n^oo 

= 1 — lim S{n) 



n— f oo 
oo 



= 1 — > r]i lim exp 

— ^ n— f oo 

1=1 

Since t{n) ~ Q{n^ 

lim exp 



{cd)^n 

we furthermore have 

Pi 



i{n) 



which gives lim„^oo P{72;(C(iY^) < t{n)} = 1. This com- 
pletes the proof for a ~ 2. 

Next, consider the case of a G (0, 2). Similarly to the proof 
for a = 2, we can prove this case by substituting r ~ Cd\pn 
and t ~ t(n) into ( fTsl l and showing that 



lim P{T,(crfy^) < £(n)} = L 



(20) 



Since the dominant decay rate |Ai| scales as 0(r^") = 
0(n^"2), by using approaches in the proof for a = 2, we 
can show (|20]|. Due to similarities, we omit the details. □ 



Lemma 4 {Lower bound for Levy flight): Suppose that un- 
der Levy flight with a distribution parameter a G (0, 2], the 
time t = t(n) in P\T(cd\fn) > t} scales as 8(rit~'^) for an 
arbitrary e > 0. Then, we have 

lim P{T(cdV^) > i{n)} = 1, 

n— s-oo 

which shows that the critical delay Cn{n) under Levy flight 
scales as ^l{n^). 

Proof: We will prove this lemma by showing that 
lim„_^oo P{Tx{cd\^/ V2) < i{n)} = 0. Then, by substituting 
r = Cd^/n and t = i{n) into ^ and taking a limit to n, we 
obtain 

lim P{T{cd^/^) < i{n)} <2 lim P{Tx{cdV^/V2) < i{n)} 

= 0. 



N denotes a set of positive integers. 



That is, lim„_^oo P{r(crf-y/n) — ^("■)} — 0' or equivalently, 

1, which proves the lemma. 



lim„^oo P{T{cdy/ri) > t{n)} 
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First, consider the case of a = 2. We substitute r ~ 
Cd\pnj \/2 and t ~ i{n) into dTsl l. Then, the series on the right- 
hand side of dTsl l becomes a function of n, and analogously 
to the proof of Lemma |3] we let 



P{T;(cdV^/V2) < i{n)} = 



tin) 



= l-S{n). 

Similarly to the proof of Lemma [3] we will show that there 
exists a constant n G N such that the infinite series S{n) 
converges uniformly on 2? = [h,oo). 

Since i{n) = <d{n^~'^), there exist constants h e N and 
c > such that 



i{ri) > cn^ 



for all n> h. 



(21) 



For a technical purpose for showing the uniform convergence, 
we restrict the domain of n as T>d = [n, d] for an arbitrary 
d> n. Let m = F{n/\/2cd)^cd~'^ . Then, the i-th function of 
the series S{n) is bounded by a constant Mi = ^{exp(— m)}* 
for all n E T>d as follows: 



Vi exp 



2p^ 



i{n) 



4 

< — exp 

TT 

4 

< — exp 

TT 

= M,. 



[cdYn 
Fi{-Kf 



cd 



Here, the first inequality comes from the bounds |77i| < ^ Vi G 
N and (l2Tl i. and the second inequality comes from the bounds 
i'^ > i Wi e N and n^" > d'" V?i £ Vd- Note that the 
series J^'iLi ^'^i converges since it is a geometric series with 
a common ratio exp(— m) e (0, 1). Hence, the infinite series 
S{n) converges uniformly on Vd- Since d is arbitrary, we get 
uniform convergence on 2?. Consequently, we can interchange 
the order of limit and summation, and we have 

lim P{T,{cdV^/V2) < i{n)} 
= 1 — lim S{n) 

2p, 



1 



OO 

E 

i=l 



^ rji lim exp ( — 

=1 

Since i{n) — <d{n^ 

lim exp 

which gives 

lim P{Tx{cd\/n 



{cd)^n 

), we furthermore have 

2p, 



i{n) 



{cdYn 



t{n) 



1, 



2) < t{n)] = 1 



4 = 1 



Note from ^ that ?{T^{cdy/^/ ^/2) < 0} = 1 - Y.T=i Vi- 
In addition, it is obvious that ?{Tx{cdy/n/ \/2) < 0} = 0. 
Therefore, we have YiuYn~foo'P{Tx{cd\/nl ^/2) < i{n)} = 0. 
This completes the proof for a ~ 2. 

Next, consider the case of a G (0,2). Similarly to the 
proof for a = 2, we can prove this case by substituting 
r ~ Cd\/n/^/2 and t = t{n) into ( fTSI l and showing that 

lim ?{Tx{cd^/n/^/2) < i{n)} = 0. (22) 



Since the dominant decay rate |Ai| scales as 6(r~") — 
Q{n~T), by using approaches in the proof for a = 2, we 
can show (|22]|. Due to similarities, we omit the details. □ 



Combining Lemmas [3] and |4] yields the following theorem. 

Theorem 1: The critical delay Cn(w) under Levy flight 
with a distribution parameter a E (0, 2] scales as 0(nt). 

Remark 1: The main idea behind the proof of Lemmas |3] 
and m was that the smallest decay rate in the exponential 
functions in (fTsl l (i.e., ^) determines the limiting behavior 
of the first exit time distribution. That is, the smallest decay 
rate characterizes the critical delay under Levy flight. 

VI. Critical Delay Analysis for Levy Walk 

In this section, we will show that the critical delay Cu{n) 
under Levy walk with a distribution parameter a scales as 
e(n^) for a £ (0, 1) and e(nf ) for a £ [1,2] (TheoremEJ. 
In Section IVI-AI we explain technical challenges and our 
approach for proving Theorem |2] In Section IVI-BI we prove 
Theorem|2]by showing that the upper bound on Cnin) scales 
as 0(n2 ) for a e (0, 1) and 0{n^) for a G [1, 2] (Lemma|6]l, 
and that the lower bound on Cnin) scales as f2(n2) for 
a G (0, 1) and n{nT) for a G [1, 2] (Lemma Hi. 

A. Technical Approach 

We first explain the technical challenges that preclude the 
use of our technique for Levy flight as well as other conven- 
tional techniques. We next explain our technical approach to 
deal with these challenges. The technical challenges are two- 
folds and are mainly inherent in the Levy walk nature. 

(i) We begin with the description of differences between 
Levy flight and Levy walk from a modeling perspective. Let 
ti{i = 1,2,...) denote the time instant when the i-th step 
begins. We take the time ti as the embedded point of the 
process {X{t)}t>o, and focus on the corresponding embedded 
process {Ei}i^fi = {X(ti)}igN- Under both Levy mobility 
models, at each embedded point ti, the destination of the next 
step of the i-th step (i.e., -Ei+i) is chosen independently of 
the past locations at time t < ti and depends only on the 
current location at time t ^ ti. That is, the embedded process 
{-EijigN satisfies the following Markov property: 

P{Ei+i = a3,+i I Ej = Xj,j 1, . . . ,i} 
= P{E,+i = Xi+i I E, = Xi}. 

Thus, under both Levy mobility models, the process 
{X{ti)}i^fi becomes a discrete-time Markov chain. However, 
the fact that the embedded point ti is chosen in a different 
way for Levy flight and Levy walk incurs the key challenge. 
In the case of Levy flight, it is chosen deterministically as 
ti = i — 1. Therefore, Levy flight is a discrete-time Markov 
process. However, in the case of the Levy walk, the embedded 
point is chosen stochastically and is correlated with step size 
as follows: ti = J2j=i (where Zj is a random variable 
denoting the j-th step size). Therefore, the Levy walk is a 
semi-Markov process [i20l whose embedded process becomes 
Levy flight. 
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Cd'fn 



Fig. 2. An example of the random variables N{n) (= 4), Zi and Z 



JV(n)- 



(ii) The proof of Lemma |2] also shows that, for a given 2- 
dimensional Levy walk, its 1 -dimensional projected processes 
also have a power-law type of step size distribution. However, 
the velocity of the projected processes is not a constant 
for every step, which implies that neither of 1 -dimensional 
projected processes of 2-dimensional Levy walk can be 1- 
dimensional Levy walk. 

Consequently, the technique used for Levy flight in this 
paper is not applicable because it requires decoupling of space 
and time. In addition, the occupation probability P{x, t) is not 
available and the derivation is not mathematically tractable. 

To cope with these technical challenges, we propose a dif- 
ferent approach based on a stochastic analysis technique char- 
acterizing the embedded Markov process of a semi-Markov 
process. Specifically, our approach is to derive a relation 
between the first exit time under Levy flight (i.e., embedded 
Markov process) and that under Levy walk (i.e., semi-Markov 
process). From this relation, our technique derives a tight 
upper bound for the critical delay. Then, by combining the 
upper bound and a lower bound for the critical delay inferred 
from our analytical result of Levy flight in Section |Vl we can 
provide the exact order of the critical delay under Levy walk. 

B. Analysis 

Let N{n) be a random variable denoting the number of 
steps occurred until t < T{cd\/n). Then, 



T{cd\/n) 



if N{n) = 1, 
E.=^^"' Z^ + ZNin) if N{n) > 2, 



(23) 



where Z^(„) is a random variable denoting the moving 
distance during the A^(7i)-th step until exiting the disc 
(See Fig. |2]). Note that ^jv(n) is not identically distributed 
with Zi, and we have 



ZN(n) < 2crf-y/n with probability I. 



(24) 



The random variable N{n) is closely related to the first exit 
time under Levy flight, denoted by Tu;{cd\/n), as follows: 



N{n) = ^TLF(cd^^)l , 



(25) 



where \x~\ denotes the smallest integer larger than or equal 
to X. In Lemma |5] we derive the order of E[N{n)], which 
will be used to study the critical delay under Levy walk. 



Lemma 5: 'E\N{n)] scales as 8(712 ) for a S (0,2]. 

Proof: From Lemma [3] we have lim„^oo 'P{Tup{cd^/n) < 
tin)} = 1 when i{n) = 9(n^+^) for a G (0,2] and an 
arbitrary e > 0. Hence, we have 



E[TLF(crfV^)] = 0(7it+^-). 



(26) 



From LemmalU we have lim„^oo ^{Tije{cd\/n) > i{n)} = 1 
when i{n) = Q{n^^'^) for a £ (0, 2] and an arbitrary e > 0. 
Thus, we have 



(27) 



By choosing e and e arbitrarily small, from ( |26] l and ( |27] l, we 
have 

E[TMcdV^)] - e(nt) Va G (0,2]. (28) 
Note from ((25) that 

E[Ti^picdV^)] < E[Nin)] < E[nF{cdV^i)] + 1, (29) 

which shows that the order of E[Af(?7)] is the same as that of 
E[Ti^F{cd\/Ti)]. Therefore, combining ( |28] ) and ( |29l ) yields the 
lemma. □ 

With the help of Lemma |5] we can derive an upper bound 
for the critical delay under Levy walk. 

Lemma 6 (Upper bound for Levy walk): Suppose that un- 
der Levy walk with a distribution parameter a, the time 
t = i{n) in ¥{T{cd^/n) > t} scales as e(n^+'^i) for an 
arbitrary ei > and a £ (0, 1), and 8(ri,'^+'2') foj- arbitrary 
£2 > and a G [1,2]. Then, we have 

lim P{T{cdV^) > i{n)} = 0, 

which shows that the critical delay Cn{n) under Levy walk 
scales as 0(71,2) for a G (0,1) and 0{n^) for a G [1,2]. 

Proof: Using Markov's inequality ||29l , we have 

E[T{cdV^)] 



P{T{cdV^) > i{n)} < 



i[n) 



(30) 



We calculate the expectation E\T{cdy/ri)\ on the right-hand 
side of ( I30I 1 by conditioning on the values of N{n) as 

E[T{cdV^i)] = E[E[^(cd^^) | N{n)]\ 

00 

= ^ E[T{cdVn) I N{n) = k] ■ ¥{N{n) 
fc=i 

From ( |23] l, we have for fc = 1, 

E[T{cdV^) I N{n) = k]= CdV^. 
In addition, from (l23l l. we have for fc = 2, 3, . . ., 

E[T{cd^/^) I N{n) = fc] 

'Af(n)-1 

^ Zi + ZN(n) Nin) ^ 
i=l 

fe-1 

= E[Z, I N{n) = fc] + E[Zk I N{n) = fc]. (33) 

i=l 



fc}. 

(31) 
(32) 
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The random variables (i = 1, . . . , fc — 1) and Zk in ( [33T l are 
correlated with the random variable N{n) {— k), whereas the 
random variables Zi {i — k + 1, k + 2, . . .) are independent of 
N{n) (= k). Specifically, for i = 1, . . . , fc — 1, the step size 
Zi should be less than the diameter of the disc 'D{n) (i.e., 
2cd\/ri). In addition, the truncated step size Zk should satisfy 
the inequality in (l24l i. Hence, the conditional expectations on 
the right-hand side of ( l33l l are bounded as follows: 



B[Z,\N{n) = fc] < E[Z\Z < 2cdV^], 
E[Zk\N{n) = k]<2cdV^, 



(34) 



where Z denotes the generic random variable for Com- 
bining ([3T|i-(l34li. we obtain an upper bound for E[T{cciV^}] 
as follows: 

E[T(cdV^)] 

oo 

< CdV^ ■ P{N{n) = 1} + 2cdV^^P{N{n) = k} 

k=2 



+ E[Z\Z < 2cdV^] ~ 1) ■ P{^W = k} 

k=2 

oo 

< 2cdVn^P{N{n) = k} 



k=l 



E[Z\Z < 2cdVn] k ■ P{N{n) = k} 



fc=i 



2cdV^ + E[Z|Z < 2cd^\ ■ E{N[n)\ 



(35) 



Using (O, we can calculate the conditional expectation 
E\Z\Z < 2cd\pr)\ in and it scales for each a G (0,2] 
as follows: 



for a e (0, 1), 
for a = 1, 
for a g (1, 2), 
for a = 2, 



E\Z\Z < 2cdM 

( a (2crfV7r)'-°~l 
1-a l-(2cdV^)-° 

\og(2cd^/ri) 
l-(2cd\Ar)-i 

a l-(2cdv/7?)^~° 
a-1 l-(2cd\AJ)-° 
^/2a cxp(-l/2g-^)~cxp(-2(cd)^n/g-^ 
, crf(cd\/2rr/o-)-crf(l/V2(T) 

e(n(i-")/2) for a G (0,1), 
6(log(n)) for a — 1, 
9(71°) for a e (1,2]. 



Since E[N{n)] scales as 8(n 2 ) by Lemma|5] the term on the 
right-hand side of (l35T l scales as 



2cdVn + E[Z\Z < 2cdV^] ■ E[N{n)] 
'e(n^) for a e (0,1), 

Q{n^ log(n)) for a = 1, 
e(n7) for a e (1,2]. 

°E[Z, I N{n) = k] = E[Z] for j = fc + 1, fc + 2, ■ ■ ■ . 



(36) 



Thus, we have from ( [35] l and ( [36] l the following: 

E[T{cdV^)] 



lim 

n— >oo 

< lim 

= 0. 



£(n) 

2cdV7I + E[Z\Z < 2cdV^] ■ E[N{n)] 
i{n) 



Therefore, from ( [30l l, we have 

lim P{r(cd0I) > i{n)} < lim ^[^(crfV^)] ^ 

n-i-oo n-i-oo t(7l) 

i.e., lim„_j.oo P{T{cd\^) > t{n)} = 0. This completes the 
proof. □ 

Lemma 7 (Lower bound for Levy walk): Suppose that un- 
der Levy walk with a distribution parameter a, the time 
t = i{n) in P{T{cd\/n) > t} scales as 0(n2~<^i) for an 
arbitrary ei > and a S (0, 1), and Q{n^~'^^) for an arbitrary 
£2 > and a G (1,2]. Then, we have 

lim P{T(cdV^) > i{n)} = 1, 

n— f 00 

which shows that the critical delay Cn{n) under Levy walk 
scales as ^{n^) for a E (0, 1) and fi(n7) for a E [1,2]. 

Proof: We will prove this lemma by showing for each of the 
cases of a e (0, 1) and a e [1,2] that 

lim P{T(crf^/H) < i{n)} = 0. 

We first consider the case of a G (0, 1). Since a Levy walker 
moves with a constant velocity u = 1, it takes at least Cd\Jn 
time to exit from the disc V^n). Thus, it is obvious that 



P{T{cd^) < CdV^i} - 0. 



(37) 



Since i{n) = 0(n2 '=1), there exists a constant n e N such 
that t{n) < Cd\pn, for n > n. Hence, we have for n > n 

P{TicdV^) < i{n)} < P{T{cdV^) < CdV^}- (38) 
Combining dJTj i and ( [38] l and then taking limits, we have 

lim P{T{cdV^i) < i{n)} < lim P{T{cdVTi) < CdV^i} 

n — ^00 n — yoo 

= 0, 

i.e., lim„_>oo P{T{cd\/n) < i{n)} = 0. We have proved the 
lemma in the case of a S (0, 1). 

We next consider the case of a G [1, 2]. In the following, 
we use the notations Tlf(-) and Tlw(-) to distinguish the first 
exit times between Levy flight and Levy walk. We will show 
based on ^ and ^ that for t > 0, 



P{Tlw(c<iV^) < t} < P{rLF(cdxAI) <t + l}. 



(39) 



From (|23ll, if N{n) = 1, then rLw(cdV^) = Cdy/n > = 
N{n) - 1. In addition, if N{n) > 2, then TLw(cdV") = 
EfJf^' Z^ + Zm> E.='r^"' > N{n)~l, where the last 
inequality comes from the assumption that the step size Z has 
a lower bound at 1 (given in Section lIVI i. Combining above 
two cases gives 

rLw(cdV^) > N{n) - 1. 
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1 2 1 2 



a a 
(a) Levy flight (b) Levy walk 

Fig. 3. Critical delays under Levy flight and Levy wallc for different a. 

From we obtain N{n) = ^^LF(cd^A^)] > Tupi^Cd^/n). 
Thus, we have with probabiUty 1, 

7Lw(cd\/n) > T]j{cdVn) - 1. 
This proves ( [39] l. Substituting t = t{n) into (|39] l. we obtain 

P{rLw(c<iV^) < t(n)} < P{^LF(c<J^/H) < + (40) 

For i{n) scaling as i{n) + 1 also scales as 

Consequently, by LemmalU the probability on the 
right-hand side of ( |40] | becomes in the limit: 

lim P{TLF(cdV^) < <(n) + 1} = 0. 

>-oo 

Therefore, from ( |40l ). we have lim„^oo P{7Lw(cd\/?i) !i 
i(7i)} = 0, which proves the lemma in the case of a e [1,2]. 
This completes the proof. □ 
Combining Lemmas |6] and [T] yields the following theorem. 

Theorem 2: The critical delay Criin) under Levy walk 
with a distribution parameter a scales as 0(n,2) for 
a 6 (0, 1) and 0(nt) for a e [1, 2]. 

VII. Discussion 

We summarize the high-level interpretations of this paper. 
Fig. [3] shows the critical delays under Levy walk and Levy 
flight, parameterized by a. Levy flight shows that the critical 
delay proportionally increases with a. However, in the case of 
the Levy walk, we can find a phase transition such that when 
a E (0,1), the critical delay is constantly 9(n2) and shifts 
to the proportional increasing phase when a G [1,2]. Two 
different scaling regions are essentially related to the fact that 
the mean step size of Levy walk for a G (0, 1) is infinite but 
finite for a G [1, 2]. In contrast to Levy walk, the travel time 
independence of step size in Levy flight leads to continuous 
scaling over a. Note that for a = 2 (i.e., Brownian motion) 
our result coincides with that in fS) which also studied the 
critical delay under Brownian motion. 

By using values of a from experimental measurements from 
Q, we can see how network delay scales with human mobility 
in practice. To give an insight to the readers, we show a 
values measured from five different sites in Table [llpresented 
in Q with a flight extraction method, "rectangle"!^ We see 

"We do not present a values from other extraction methods in (7) which 
intentionally exclude some detailed motions of real traces. To capture specific 
behaviors of humans, one can borrow those a values. 



TABLE I 

Experimental Q values for different sites presented in Q. 



Site 


a 


Site 


a 


KAIST 


0.53 


New York City 


1.62 


NCSU 


1.27 


Disney World 


1.20 






State fair 


1.81 



that critical delays for human mobility range from 0(n°-^^) to 
Q{nP-^^). Human mobility mainly has a > 1, in which case a 
longer delay than O(y^) is needed. This implies that it may 
be hard to design a low delay protocol for mobile networks 
under human mobility. 

Our contribution is not restricted to the mathematical deriva- 
tion of delay scaling for new mobility models. We provided 
techniques that connect the diffusion equation of a continuous 
time random walk to the delay scaling as well as that analyze 
the delay scaling of semi-Markovian movements. We expect 
that our techniques can be further developed to the analysis 
of other detailed performance metrics such as contact time 
distribution and the generalized delay-capacity tradeoff for 
various levels of per-node throughput. 

Future work includes investigation of throughput and delay 
scaling for mobile networks with heterogeneous and collective 
node mobilities. In addition to the recent research topics on 
"per-node throughput scaling" under inhomogeneous spatial 
node distributions (i.e.. Cox process, Neyman-Scott process, 
Matern cluster process and Thomas process), e.g., 1301 , ISTI , 
our paper can be an important step to the study of delay 
scaling under such heterogeneous networks. There is an insight 
from [8 1 that in human-assisted networks, the actual delays 
might be even shorter This is because human mobility is not 
completely random: people tend to visit the same locations 
and regularly meet a group of people every day. Although their 
mobility can be characterized by heavy-tail distributions, these 
regularity in daily mobility significantly facilitates routing of 
packets among people (as long as they are socially connected). 
Therefore, there remains a possibility of designing a low delay 
protocol for mobile networks under heterogeneous human 
mobility by judiciously utilizing these social factors. 

VIII. Conclusion 

We have presented Levy mobility models consisting of Levy 
flight and Levy walk parameterized by a and studied the 
critical delay under both mobility models. Levy mobility is 
known as a realistic human mobility so that the critical delay 
we provided here can be essential in designing an architecture 
and protocols of a wireless mobile network. The insight that 
the critical delay scales as Q{n^) for Levy mobility models 
in the range of a e [1,2] is especially important because 
it is anticipating that the delay of mobile networks with 
human mobility (e.g., smartphone networks, pocket switched 
networks) could be quite high in practice, considering the a 
values measured in real traces. The insight tells that mobile 
networks operated by human mobility patterns may need to 
prepare an alternative path for delay sensitive data as well as 
even for delay tolerable data whose tolerance level is limited. 
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Appendix A 

Critical Delay Analysis for Levy Flight without 
Continuous Limit 

In Section |V] we have studied the critical delay under 
Levy flight using continuous limit. By following the technique 
in 161, we can study the critical delay without continuous 
Umit (i.e., with a non-zero scale factor |c|) and can derive a 
lower bound for the critical delay under Levy flight. Lemma [H] 
summarizes the result. 

Lemma 8: With a non-zero scale factor |c|, the critical delay 
Cn(n) under Levy flight with a distribution parameter a G 
(0,2] scales as il{n'^). 

Before proving the lemma, we give a remark. The scaling 
property of the critical delay with continuous limit (shown 
in Theorem [T]) works as an upper bound for the one without 
continuous limit. Hence, the result in Lemma [8] shows that our 
analysis in Section [V] gives the tightest upper bound, which 
justifies our technique using continuous limit. We now give 
the proof of Lemma [8] 

Proof: Similarly to the proof of Lemma ID we will prove this 
lemma by showing that 

lim P{r,(cdV^/V2) < i{n)} = 0, (41) 

where t{n) = for an arbitrary e > and a e (0, 2]. 

Then, from we obtain lini„_j.oo P{7'(cd\/?i) — K^)} — 
21im„^ooP{^I:(cdV"/^/2) < = 0. That is, we 

have lmln-^co P{T{cdy/n) < t{n)} = 0, or equivalently, 
linln-^ocP{T (cd^/n) > t{n)} = 1, which shows that the 
critical delay Cuin) scales as 51(ri'5). 

Without loss of generality, we assume X(0) = (0, 0). Then, 
from (|9]l, Xx{t) for t = 1,2, . . . can be expressed as 



x,{t) = Y,z, cos e,. 

1=1 



(42) 



Let Zcos9 denote the generic random variable for Z, cos0j. 
By the independence of random variables Z and 6, the 
mean of Zcosd is given by E[Zcos9] ~ E[Z]E[cos9] = 
0. Hence, from (gUi, the mean of Xx{t) is given by 
E[X2;(<)] = tE[Zcos9] = and the variance of Xx{t) 
becomes E[(X^(i))2] = tE[{Zcose)% Thus, Hoeffding's 
inequality |j6l gives an upper bound for P{X4t) > r/y/2} 

as follows: P{Xx{t) > r/V2} < cxp (- 8,e1(/cos.)^] ) • By 
the symmetry of node motion, we have 

P{\Xx{t)\ > r/V2} < 2cxp f- ,,^r/ ^^^91 ^ ■ (43) 



^ 8tE[{Zcosey]y 

Due to (Al), the event {T^(r/^/2) < k} for k = 1,2,.. 
implies the event Ut=i{l^a:(OI — ^/V^}. Hence, we have 

fe 

P{T4r/V2) <k}< ^P{|X,(i)| > r/V2} 
t=i 



< 2 ^ cxp - 

t=i 



< 2k exp 



8tE[(Zcos6')2 

^2 



8fcE[(Zcos( 



, (44) 



where the second inequality comes from ( |43] ). Substituting r 
Cdy/n and k = t{n) into (l44l i. we have 



P{Tx{cdV^/V2) < t{n)} 
< 2i{n) exp ( ^ 



{cdfn 



8t(n)E[(Z cos 61)2] 



(45) 



In the following, we will derive a bound for E[(Zcos0)^]. 
Since E[(Zcos6l)2] = E[(Z| cos 6*1)2], we have 

E[(Zcos0)2] = / z''AFz\,osep)- (46) 
Jo 

We first consider the case of a G (0,2). From the CCDF of 
Z\ cos 6*1 given for z > 1 in (fT2b . we have for z > 1, 



■ ■ — = P{Z\ cosy > z\ 

dz dz 



< 



2ac{n) 
ac*c{n) 



(cosi?)"di5 



< 



Thus, the integral on the right-hand side of (|46] | is bounded 
above by 

z'^dFz\cose\{z) 

r-l 

= J dFz\ cos 0\{z) + J z'^dFz\cose\iz) 
P{0 < Z\ cos9\ < 1} + J^^ z^^^^dz 

= P{0 < Z\ cos9\ < 1} + °'g*^(") („i-f _ 

2 — a 

from which we have 

E[(Zcos6')2] = 0(ni~t) for a e (0,2). (47) 

We next consider the case of a = 2. By following the approach 
in the case of a € (0, 2), we have 

E[(Zcos6l)2] = 0(7ii-t) for a = 2. (48) 

Combining (gT) and (gUl gives E[(Z cos 61)2] ©(n^-f) for 
a G (0, 2]. Hence, there exist constants n G N and c > such 
that 

E[(Z cos 61)2] < cn^-t for all n > fi. (49) 

In addition, since t{n) = 8(ri,'^~'^), there exist constants h E 
N and c > such that 



i{n) < cn'^ ^ for all n> h. 



(50) 



By ( |49] l and ( |50] l, the term on the right-hand side of ( |45] ) is 
further bounded by 

^ ^ 8t(n)E[(Zcos0)2], ^^^^ 

< 2cn 2 exp ' 
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By L'Hospital's rale, ( fSTT l becomes in the limit as 



lim 2t{n) cxp 



< lim 2cn ^ exp 



8t(n)E[(Z cos 61)2 



(52) 



0. 



Combining ( |45] l and ( |52] i proves ( |4TI ). This completes the 
proof. □ 
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